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Weakly-coupled systems in quantum control 

Nabile Boussaid, Marco Caponigro, and Thomas Chambrion 



Abstract 

Weakly-coupled systems are a class of infinite dimensional conservative bilinear control systems 
with discrete spectrum. A property of these systems is that they can be precisely approached by finite 
dimensional Galerkin approximations. This feature is of particular interest for the approximation of 
quantum system dynamics and the control of the bilinear Schrodinger equation. 

This paper provides rigorous definitions and analysis of the dynamics of weakly-coupled systems 
and gives sufficient conditions for an infinite dimensional quantum control system to be weakly-coupled. 
As an illustration we provide examples chosen among common physical systems. 

Index Terms 

Quantum system, Schrodinger equation, bilinear control, approximate controllability, Galerkin ap- 
proximation. 

I. Introduction 

A. Physical context 

The state of a quantum system evolving on a finite dimensional Riemannian manifold f2, with 
associated measure //, is described by its wave function, that is, an element of the unit sphere 
of L 2 (Q, C). A system with wave function ip is in a subset oj of Q with probability \tjj\ 2 diJ,. 
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When a system is submitted to excitations by p external fields (e.g. lasers) the Schrodinger 
equation reads 

Bib 1 p 

where A is the Laplace-Beltrami operator on f2, V : — > R is a real function, usually called 
potential, carrying the physical properties of the uncontrolled system, Wi : Vt — > R, 1 < / < p, 
is a real function modeling a laser /, and ui,l < I < p, usually called control, is a real function 
of the time representing the intensity of the laser I. 

A natural question, with many practical implications, is whether one can find a control 
(wi, • • • , Up) for which the associated dynamics generated by © has the desired behaviour (for 
instance, steering some given initial data to some given target). In the case p — 1, considerable 
efforts have been made to study the controllability. We refer to|[D,ll3,0.lffl.l0.Gl.[D.lllD. 
ED, iflOl and references therein for a description of theoretical results concerning the existence 
of controls steering a given source to a given target or a neighborhood of it. 

B. Finite dimensional approximations 

The main difficulty in the study of (OQ) is the fact that the natural state space, namely L 2 {Vt, C), 
has infinite dimension. To avoid difficulties when dealing with infinite dimensional systems, for 
example when studying practical computations or simulations, one can project system CD on finite 
dimensional subspaces of L 2 (f2, C). A vast literature is now available on control of bilinear finite 
dimensional quantum system, we refer for instance to [fTT]|. |fT2|. A crucial issue is to guarantee 
that the finite dimensional approximations have dynamics close to the one of the original infinite 
dimensional system. 

As a matter of fact, a special class of bilinear systems of the type of ([I]) which are said 
weakly-coupled (see definition in Section UU) exhibits very nice properties of approximations (see 
Proposition H] below). The notion of weakly-coupled system has been used in |[T3l implicitly, in 
the case p — 1 and W bounded. 

The aim of this work is to provide an analysis of weakly-coupled systems, to present a 
sufficient condition for controllability for these systems, and to show that two important types 
of bilinear quantum systems frequently encountered in the literature are weakly-coupled. 
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C. Content of the paper 

In Section [II] we define weakly-coupled systems, we state some properties of their finite 
dimensional approximations, and we prove a controllability result. We then study two important 
examples of weakly-coupled systems, the first (Section HITl) covering, among others, the case 
where Vt is compact and the second (Section HVT) the case where the system is tri-diagonal. 

II. Weakly-coupled systems 

A. Abstract framework 

We reformulate CQ) in a more abstract framework. This will allow us to treat examples slightly 
more general than (OQ), for instance, the example in Section IIII-Al In a separable Hilbert space 
H endowed with norm || • || and Hilbert product (•,•), we consider the evolution problem 

^ = (A + X>W (2) 
where (A, Bi, . . . , B p ) satisfies Assumption [TJ 

Assumption 1. (A, B x , . . . , B p ) is a (p + l)-uple of linear operators such that 

1) for every u in R p , A + Y^i^iBi is essentially skew-adjoint on D(A) and i(A + Y^iUiBi) 
is bounded from below; 

2) A is skew-adjoint and has purely discrete spectrum (— iX^keN, the sequence (\ k )keN is 
positive non-decreasing and accumulates at +oo. 

In the rest of our study, we denote by (4>k)keTsi an Hilbert basis of H such that A<p k = —\\ k (f) k 
for every k in N. We denote by D(A + J2i uiBi) the domain where A + J2i u iBi is skew-adjoint. 

Assumption 1 1 1 1 1 ensures that, for every constants ui,...,u p in R, A + J2i u iBi generates a 
group of unitary propagators. Hence, for every initial condition ^ m H, for every piecewise 
constant control u : t G R — > J2n=o un X[t n ,t n +i) (*) e ~& p wu " n — *o < *i < • • • < t N+1 we can 
define the solution of © by 1 1-> T^ijj , where 

T« = e (t-t 3 -_i)(^+E«f" lB o e (* J -i-< J -2)(^+E«(" 2B i) o • • • o e *o(^+E«? B i) ; 
forte fe-i, tj), j — 1, . . . , N. 
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Remark 1. From Assumption 1111 1 we deduce that for every u G W, A + Y,i u iBi is bounded 
from D(A) to H as well as J2i UiBi. As a consequence, the resolvent of A + J2i u iBi is compact. 
Since i(A + J2i u iBi) is bounded from below, the spectrum accumulates at +oo. 

B. Energy growth 

From Assumption II 121 the operator iA is self-adjoint with positive eigenvalues. For every -0 
in D(A), \Aip = EjeN ^j(4>ji i>)$y F° r every s > 0, we define the linear operator \A\ S := (iA) s 
by \A\ s ij = E ieN Xj{M)<l>3, for every V in £>(|A| S ) — {ip G H : £ ieN Af |(<^>| 2 < +oo}. 
We define the s-norm by \\ip\\ s = |||^4| s, 0|| f° r every tp in Z)(|t4| s ). The 1/2-norm plays an 
important role in physics; for every ip in Dd^l 1 / 2 ), the quantity \(Aip,ip)\ = \\ip\\l/ 2 i s me 
expected value of the energy. 

The notion of weakly-coupled systems is closely related to the growth of the expected value 
of the energy. 

Definition 1. Let He a positive number and let (A, Si, ... , B p ) satisfy Assumption 11111 Then 

(A, B h . . . , B p ) is k-weakly-coupled if for every u G R p , D(\A + EiUiB^/ 2 ) = D(\A\ k ' 2 ) 
and there exists a constant C such that, for every 1 < I < p, for every ip in D(\A\ k ), 
m\A\ k ^,Brt)\<C\(\A\ k iP^)\. 
The coupling constant Ck(A, B\, ... , B p ) of system (A, Bi, ... , of order k is the quantity 

\m\A^j^m 
SU P SU P l/l /life/ /\l ■ 

We have the following technical interpolation result whose proof is given in the appendix. 

Lemma 1. Let A and A' be invertible skew-adjoint operators with compact resolvent. Let k be a 
positive real Assume that D(\A\ k ) = D(\A'\ k ). Then for any real s G (0, k), D(\A\ S ) = D(\A'\ S ). 

A first property of the propagator of a weakly-coupled system is given by the following 
proposition. 

Proposition 2. Let k be a positive number and let (A, Bi, . . . , B p ) satisfy AssumptionUland be k- 
weakly-coupled. Then, for every ip G D(\A\ k / 2 ), K > 0, T > 0, and piecewise constant function 
u = (u u ...,u p ) for which £f =1 < K, one has ||T^(^ )|| fc / 2 < e Cfc ( A ' Bl >-'^ x ||^ ||fc/2- 
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Proof: We present here a simple proof in the special case where _D(|/l| fc+1 ) = D(\A + 
J2i< p uiBi\ k+1 ) for every u in R p . This equality holds for the most common physical examples. 
A general proof of Proposition [21 involving rather advanced regularization techniques to relax 
this extra assumption is presented in the Appendix. 

First note that, for every t > and for every (ui, . . . ,u p ) in R p , the set /}(|v4| fe+1 ) = 
D(\A + J2iUiBi\ k+1 ) is invariant for the unitary map ip i-)- e t<yA+ ^i UlBl ^ip. Moreover, for every ip 
in D{\A+Ei uiBi\ k+1 ), the map \A+Ei uiBi\ k e t{A+ ^i UlBl) ip = e' (j4+ S; u i B i)\A+J2i uiBi\ k ip 
is C 1 from [0, +00) to H, with derivative t m- (A + J2i u l B l )e t{A+ 'E' UlBl) \A + EiUiB t \ k ip = 
\A + EiUiB t \ k (A + J2iUiBi)e t{A+ ^' UlBl) ^j. In other words, the map t \-> e *< A+ £» ttjB, >V is C 1 
from [0, +00) to D(\A + E/ «i-B,| fc ) = Z?(|^4[ fe ). 

Fix w : [0, +00) — > R p piecewise constant, ip in D(|y4| fe+1 ) and consider the real map 
/ : t h> (|A| fc T^o, T^Vo)- Since ip belongs to £>(|A + Ef=i then / is absolutely 

continuous and for the argument above is piecewise C 1 . For almost every t, 

= <|A|*T?^o, + E»l(*)B,)TrV'o> + <|>t|*(j» +E«,(<)Bi)T?V'o,TrVo) 

1=1 (=1 

= 2K(|A| fc T^o, (^4 + E Wz(t)A)T^o) 

z=i 

1=1 

Since (A, Si, ... , B p ) is /c-weakly-coupled, then 

|/'(*)| < ^J2\ui(t)\\(\A\ k T^ ,BiT^ } \ < 2c k (A,B ll ...,B p )J2\ui(t)\f(t). 
1=1 1=1 

From Gronwall's lemma, we get (\A\ k ip(t),ip(t)) < e 2ck(A '- Bl '-'^ ,) ^=i^ |l " l(T) *"||V'o||fc/i! for 
every tp in D{\A\ k+l ). 

k 

Let (^o)neN be a sequence in D(\A\ k+1 ) converging to ip in D(\A\ 2) for the fc/2-norm. For 
every ip in D(|A|2j and for every t > 0, 

K|A|i^T?«)>| < k^|a|4t»(^)>i < ||^||||T?«)||j, 

which is bounded uniformly with respect to n from the first part of the proof. Since ip !->■ T^(ip) 
is unitary, (T"( , 0o))neN converges to T"(^ ) (for the norm of iJ) and \(\A\ 27/j, T"^)) I < 
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e 2c k (A,B 1 ,...,B,)jf||^|| ||^, o || Hence? T^(^o) belongs to D(\A\%)* = D(\A\%) and ||T^ )||| < 
e c k {A,B lt ...,B p )K^ o ^ k ^ which concludes the proof in the case in which D(\A\ k+1 ) = D(\A + 



Y.i< p niB l \ k + l ). m 

C. Good Galerkin approximation 

In this section we show that a weakly-coupled system admits a finite dimensional approxima- 
tion whose trajectories are close, at every time, to the solutions of the original infinite dimensional 
system. For every iV in N, we define the orthogonal projection 

n N : i/j e H i-)> E (<A?> e H. 

j<N 

Lemma 3. Let k be a positive number, (A, B\, ... , B p ) satisfy Assumption UI and be k-weakly- 
coupled. Assume that there exist d > 0, < r < k such that < d||^|| r / 2 for every ip in 

D(\A\ r ' 2 ) and I in {1, . . . ,p}. Then, for every K > 0, n E N, iV e N, (V>j)i<i<n in D(\A\ k ' 2 ) n , 
and for every piecewise constant function u = (ui, . . . , u p ), 

E IKHzi < K =► ^(Id-Tr^T^)!! < rfAS;f )/2 e Cfc ( A ' Bl ""^^||^|| fc/2 , (3) 

m=l 

/or every £ > 0, £ = 1, . . . , p and j = 1, . . . , n. 



Proof: Fix j £ {1, . . . , n}. For every N > 1, one has 

oo 

||(Id-^)T?(^)||^= E A;|(0 n ,T^,))r<A^ill^)ll' /2 - (4) 

n=N+l 

By Proposition El ||T^(^-)||L 2 < e 2ck ^ Bl '- Bp)K \\i)A\l /2 for every t > and it of L^norm 
smaller than fC Equation © follows as, for every I — 1, ... ,p, \\Biip\\ < d\\\A\^ip\\. ■ 

Remark 2. Since r < k, then ||-B;(Id — 7rjv)T"(^j)|| r ^ 2 tends to 0, uniformly with respect to u, 
as N tends to infinity. 

Definition 2. Let iV G N. The Galerkin approximation of © of order iV is the system in iJ 

2=1 

where A^ = ttnAttn and = ttnBittn are the compressions of A and 5; (respectively). 



We denote by XV N Jt, s) the propagator of ([£jv]) for a p-uple of piecewise constant functions 

U = (Ui, . . .,u p ). 
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Remark 3. The operators and b[ n ^ are defined on the infinite dimensional space H. 

However, they have finite rank and the dynamics of (Sat) leaves invariant the iV-dimensional 
space Cn = span 1<KA ,{<^}. Thus, (Sat) can be seen as a finite dimensional bilinear dynamical 
system in C N . 

Proposition 4 (Good Galerkin Approximation). Let k and s be non-negative numbers with 
< s < k. Let (A, Bx, ... , B p ) satisfy Assumption [7J and be k-weakly -coupled. Assume that 
there exist d > 0, < r < k such that \\Biip\\ < d\\ip\\ r /2 for every ip in D[\A\ r l 2 ) and I in 
{l,...,p}. Then for every e > 0, K > 0, n E N, and {^^i-cj-cn in D(\A\ k ^ 2 ) n there exists 

]VgN such that for every piecewise constant function u = (%, . . . , u p ) 

v 

J2 IMIn < K WT^tpj) - X^ N) (t, 0)71^1^2 < e, 
i=i 

for every t > and j = 1, . . . , n. 

Proof: Consider the case s = 0. Fix j in {1, . . . , n} and consider the map t >->■ n N T^(^j) 
that is absolutely continuous and satisfies, for almost every t > 0, 

ar 1=1 i=l 

Hence, by variation of constants, for every t > 0, 

ttatT^^) = Xfo(t, 0)7^ + £ /xfofc ^vr^Id - -K^T^u^dr. (5) 

i=i 1/0 

By LemmaH the norm of t ^ ^(Id - 7r JV )TJ t (^) is less than rfA^^e^^^'-'^^H^Hfe^. 
Since X^Jt, s) is unitary, 

IKatT^^) -Xfafi, 0)^-11 < ^Afcf^e^^'-^^ll^-llfc/a. (6) 

Then 

\\T^j) -Xfa(t, 0)^-11 < IKld-Tr^T^OH + IkivTr(^) - 0)^11 

(7) 

This completes the proof for s = since \ N tends to infinity as N goes to infinity. 

Note that, if X is a set and (M n ) ne N is a sequence of functions from X to H that tends 
uniformly to (the null function) for the s-norm and is uniformly bounded for the /c-norm for 
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s < k, then (w n ) nG N tends uniformly to in the ^rp-norm. This is a consequence of Cauchy- 
Schwarz inequality, indeed 

||«n||l+fc = \(\A\^U n , \A\^~U n )\ = \(\A\ S U n , \A\ k U n )\ < \\u n \\ s \\u n \\ k . 



We conclude the proof in the general case s > applying this interpolation result, combined 
with a bootstrap argument, on the sequence (un) /vgn with un '■ (t,u) i-> (Xfo {t, Q)n N — T^)ipj, 
defined on X = [0, +oo) x {u G L l : \\u\\ L i < K}. ■ 

Remark 4. In the case s = 0, there is an explicit estimate for the order of the Galerkin 
approximation which existence is stated in Proposition |4] For instance, by ©, \\Ti N T^{ipj) — 
Xy , N s(t,0)TTN'ipj\\ < £ if N is such that 

( Kde^ A ^-^ K U 3 \\ k/2 \ & 

A/V+l > I I , (P) 

for j = l,...,n. 

D. Approximate controllability in s-norm 

The a priori bounds on the A;-norm for the solution of a system are a deep obstruction to 
exact controllability, but provide powerful tools for the study of the approximate controllability. 

Definition 3. Let {A, B) satisfy Assumption [Q A subset S of N 2 couples two levels j, k in N, 
if there exists a finite sequence y(s{, si), . . . , (sf, s|)^ in S such that 

(i) s\ = j and s| = fc; 

(ii) s 2 = s{ +1 for every 1 < j < q — 1; 
{Hi) (<f> s{ ,B<f> 4 )^0forl<j<q. 

The subset S 1 is called a connectedness chain for (A, £>) if 5 couples every pair of levels in N. 
A connectedness chain is said to be non-resonant if for every (s 1; s 2 ) in S 1 , | A S1 — A S2 1 ^ | A tl — A 4a | 
for every (ti,t 2 ) in N 2 \ {(si,s 2 ), (s 2 ,si)} such that (<f) t2 , B(f) tl ) ^ 0. 

Definition 4. Let (A,£>) satisfy Assumption Q] and s > 0. The system {A, B) is approximately 
simultaneously controllable for the s-norm if for every T G U {H) (unitary operators acting on 
H) leaving invariant, tpx, . . . , tp n G Z^(|^4| s ), and e > 0, there exists a piecewise constant 

function u £ : [0,T e ] — >■ R such that 

||f^-T£^||.<e. 
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for every j = 1, . . . , n. 

Proposition 5. Let k be a positive number. Let (A, B) satisfy Assumption^ be k-weakly-coupled, 
and admit a non-resonant chain of connectedness. Assume that there exist d > 0, < r < k such 
that \\Bip\\ < d\\ \A\%ip\\, for every ip in D(\A\%). Then (A, B) is approximately simultaneously 
controllable for the norm \\ ■ \\ s /2 far every s < k. 

Proof: Fix e > 0, ip u . . . ,ip n e D(\A\ S > 2 ), and f G U(H) such that f (^i), . . . , f(ip n ) G 
D{\A\ a ' 2 ). Fix m sufficiently large such that ||f (ipA - vr ni f (^)|| s/2 < e/3 for every j = 
l,...,n. 

There exist li, . . . ,l n such that £ i-> (e liA 'i , . . . , e ttx ' n ) is e-dense in the torus T n (see ifTOl 
Proposition 6.1]). Call m = max{ni, h,..., l n }. 

By IfTOl Remark 5.9] there exists K\ such that for every 77 > there exist a control u[ satisfying 
IKi IU 1 < Ki and 9 1 , 6 2 , ■ . . , 9 n , such that ||T T * (ipj) — e 1 ^^. \\ < 77, for every j — 1, . . . , n. 

Similarly, since the hypotheses of IfTOl Remark 5.9] apply to the system (—A, —B) (see IfTOl 
Section 6.1]), we have the existence of K 2 such that for every 77 > there exists u 2 satisfying 
11^2 1 1 l 1 < K2 and 0~i, . . . , 9 n 6 R such that the solution of the system 

^(t) = -(A + u(t)B)m 

at time T 2 with initial condition T(^) and corresponding to the control u 2 is ?7-close in the 
norm of H to e 10J '</>/ . for every j = 1, . . . , n. 

Let r such that \\e TX ^e^ - e w ^ \\ < rj for every j = 1, . . . , n. Let T = T x + r + T 2 and let 
u : [0, T] — > R be the piecewise constant control denned by 

' u\{t) te^T,), 
u\t) = < £e[T 1 ,T 1 + r), 

^(T 2 -(t-T!-r)) te [Tx+T,T], 

The control vP above steers a solution of tp = (A + uB)ip with initial condition ipj 377-close in 
the norm || ■ || to Y(tpj) in a time T. 

Let K = K\ + K 2 . By Lemma [3l we have that there exists N = N(e, K,s) > n such that 

\\U\\ L 1<K =► ||(Id-7T Ar )T«(^)|| S /2<|, 

for every j — 1, . . . , n and £ > 0. 
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Note that, on span{0 1; . . . , 4>n}, we have || ■ || s / 2 < A^ 2 || • ||. Therefore for every j = 1, . . . , n, 

|f (Vi) - (^)|| s / 2 < ll(w - 7r N )(t (^.) - r^))L/2 + lkiv(f (^-) - T^))l| S /2 

< || (M - tt^T (^)|| s/2 + ||(Id - ir N m\4>j))\U/2 + 4 2 ||t(^-) - 



< j + 3A^ r] < e, 



for r] sufficiently small. 



III. The bounded case 



k 

Proposition 6. Let k be a positive integer. Assume that for every u G R p , Z)(|A|a) = D(\A + 
J2i u iBi\z) and that for every I = 1, . . . ,p the restriction of Bi to D(\A\z) is bounded for the 
—norm. Then (A, Bi, . . . , B p ) is k-weakly-coupled. 

Proof: For every I = l,...,p, let ||S^|U/ 2 < Cz,fc|MU/2 for every ip in D(\A\ k ). Then 
\(A k i),B^)\ = \(\A\ty,\A\$Bri>)\ < |||A|mill|A|^ll < C ltk \\\A\^f = d, k \(A^^)\ 
for every ip in _D(|A| fc ). ■ 

A. Example: single trapped ion 

This example is a model of a single ion with two possible states {excited state and ground 
state) submitted to a superposition of external fields. It has been extensively studied (see for 
example OH, 031, and ifToll). 

The state of the system is (ij) e ,ip g ) in H = L 2 (H, C) x L 2 (R, C). The dynamics is given by 

i^f- = u(-A + x 2 )ip e + Qip e + (u^t) cos(V2r)x) + u 2 (t)sm(V2r]x)) i[) g 

i^f = lo(-A + x 2 )ip g + n-ipg + ( Ul (t) cos(V2rjx) + u 2 (t) sm(V2rjx)) 
where u,Q,r) are positive constants related to the physical properties of the system. The two 
real valued controls u\ and u 2 are usually a sum of periodic functions with positive frequencies 
Vt, Vt + to and — to. With our notations, the dynamics reads 

d ^ = A,p + Ul (t)B^ + u 2 (t)B 2 (iP) (9) 

where A is the diagonal operator A : (ip e , ijj g ) !->■ — i(o>(— A J t-x 2 )ip e +£l?p e , to{— A+x 2 )?/> g + f2^ 9 ), 

B x : (if) e ,i>g) -i(cos(y/2rjx)ij g ,cos(V2r]x)ij e ), and B 2 : (il> e ,i> g ) ^ -i(sm(y/2r]x)^ g , sm(^/2r]x)i/; e ). 
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By IfTTl Theorem XIII.69 and Theorem XIII.70], the operator A is skew-adjoint with discrete 
spectrum and admits a family of eigenfunctions which forms an orthonormal basis of H. Since B x 
and B 2 are bounded then, for every real constants U\ and u 2 , A+uiB 1 +u 2 B 2 is skew-adjoint with 
the same domain of A (see [fT8l Theorem X.12]). The spectrum of A is the sequence (— iA n ) n£ N = 
—i(u(n+l/2)+Q) n g^. For every n in N, the eigenvalue — iA n has multiplicity 2 and is associated 
with the 2-dimensional subspace of L 2 (R, C) x L 2 (R, C) spanned by {(/ n ,0), (0,/ n )} where 
/„ is the n th Hermite function. Assumption Q] is then verified. Since, for every k in N, all 
derivatives up to order k of x >->■ cos(\/2^) and x i— > sm(y/2r]x) are bounded for the L 00 -norm 

k k 

by Ck = 22r] on R then B\ and B 2 are bounded by 2 Cj on for every k. Moreover 

for every (u\,u 2 ) G R 2 , D(A h ) = D((A + U\B\ + u 2 B 2 ) k ). Indeed by induction on k 

D{(A + u l B 1 + u 2 B 2 ) k+1 ) = G + uiSi + M 2 5 2 ) fe ) : 

(A + Ul B x + m 2 5 2 )^ G + u 1 B 1 + u 2 B 2 ) k )} 

= {iPe D(A k ) : (A + u l B 1 + u 2 B 2 )iP G ,D(A fc )} = D(A k+1 ), 

since (mBx+^B^ip G when ^ G J D(A fc ). Lemma Q] provides D(|A| S ) = L'(|A+mi j B 1 + 

n 2 _B 2 | s ) for any s > 0. Hence, by Proposition [6] the system (A, Bi, B 2 ) is /c-weakly-coupled for 
every k, with coupling constant smaller than 2 k C^. 

B. The case of a compact manifold 

We focus here on the case where the space Vt is a compact Riemannian manifold without 
boundary. By Rellich-Kondrakov and Weyl theorems, if V is essentially bounded the operator 
A = — i(A + V) : H 2 {VL) L 2 (f2, C) has purely discrete spectrum (— iA n ) n€ N with A„ non- 
decreasing to infinity (see for instance |fT9l Theorem 7.2.6]). Note that Ai is not necessarily 
positive but this is the case considering A + i{\\ — 1) instead of A. This shift gives a physically 
irrelevant phase term, ^(Ai- 1 ^ on the dynamics associated with A. 

Lemma 7. Let k be a positive integer, Q be a compact Riemannian manifold, V : Q — > R be 
C 2k {9). Then the domain of the operator (A + V) k is H 2k (Vt). 

Proof: Since Vt is compact it is sufficient to prove the proposition on a bounded domain 
of R™. The operator — L4 = A + V is an elliptic operator of order 2. By [|20l Theorem 8.10] if 
Af G H k (tt) then / G H k+2 (VL) and by induction we have that D(\A\ k ) = H 2k (VL). ■ 
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Proposition 8. Let k be a positive integer, Vibe a compact Riemannian manifold, V, W : Q — > R 
be two C 2k (Q, R) functions on Q. Define A = -i(A + V) : D(A) — >■ L 2 (fi, C) and 5 = iW : 
L 2 (f2, C) — >■ L 2 (Q, C). 77?en (A -8) w k-weakly-coupled. 

Proof: Note that for every / G C 2fe there exists a constant C fc = 2 2k+1 snp 0<j<2k \\W® \\L™(n.n) 
such that || W^/lln 2 * < Cjfc||/||ii" 2fc - From Lemma U\ the norm || ■ ||jj-a* and the A;-norm are 
equivalent. Therefore, by Proposition [6l the system is A;-weakly-coupled. ■ 

Remark 5. As consequence of Lemma [7] and Proposition [8] we have that, in the case of a compact 
manifold, if the potentials are in C m (VL) then Proposition |4] applies with k = m/2 — 1 and r = 0. 

C. Example: orientation of a rotating molecule in the plane 

We consider a rigid bipolar molecule rotating in a plane. Its only degree of freedom is the 
rotation around its centre of mass. The molecule is submitted to an electric field of constant 
direction with variable intensity u. The orientation of the molecule is an angle in Vt = SO (2) ~ 
R/27rZ. The dynamics is governed by the Schrodinger equation 



dt 

Note that the parity (if any) of the wave function is preserved by the above equation. We 
consider then the Hilbert space H — {ip e L 2 (f2, C) : ip odd }, endowed with the Hilbert 
product (/, g) = f n fg. The eigenvalue of the skew-adjoint operator A = ij^ associated with 
the eigenfunction <pk ■ 9 >-» sm(k9)/y/ir is — iA& = — iA; 2 , k e N. The domain of \A\ k is the 
Hilbert space H k = {ip e H 2k ({l,C) : ip °dd }. The skew-symmetric operator B = — icos# 
is bounded on D{\A\ k l 2 ) for every k. By Proposition [6l for every k in N, (A,B) is fc-weakly- 
coupled. Proposition |4] applies for every k with r = and d — 1. In Section HV-CI we also give 
an estimate on the coupling constant Ck(A, B) for this system. 

From the point of view of the controllability problem, notice that the operator B couples only 
adjacent eigenstates, that is (4>i,B(fij) = if and only if \l — j\ > 1. Since A; + i — A/ = 2Z + 1 
then {{j, I) e N 2 : |/ — j\ = 1} is a non-resonant connectedness chain for (A, B). Therefore, 
by Proposition \5\ the system provides an example of approximately simultaneously controllable 
system in norm H k (£l) for every k. Note that, since the eigenstates belong to H k (Q) for every 
k then the reachable set from any eigenstate is contained in H k ({l) for every k. 
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D. Example: orientation of a rotating molecule in the space 

We present the physical example of a rotating rigid bipolar molecule. Unlike last example 
the motion of the molecule is not confined to a plane. The model then can be represented by 
the Schodinger equation on the sphere. In this case, Vt = S 2 is the unit sphere, the family 
(Y^ m )^>o,| m |<^ of the spherical harmonics is an Hilbert basis of H = L 2 (tt, C), and the control is 
represented by three piecewise constant functions u\, u 2 , u 3 . The controlled Schrodinger equation 
is 



i 7^ — = (—A + Ui(t) cos 9 sin v + u 2 (t) sin 9 sin v + u 3 (t) cos v)if)(v, 9, t), (is, 9) G S 2 . 

Therefore, since Vt is compact, Proposition @] applies for every integer k with d = 1 and r = 0. 



We deal with the case where p — 1 and B couples only adjacent levels of A. 
A. Tri-diagonal systems 

Definition 5. A system (A, B) satisfying Assumption [TJ is tri-diagonal if for every j, k in N, 
| j — A; | > 1 implies ((f) j, B(p k ) = 0. 

In the following, we denote b j:k = ((f) j, B(f> k ). 



R, D((A + uB) k ) = D(A k ). Moreover, D(A k ) is invariant for e*^ 4 "" 5 ) for any u in R and t 
in R 

Proof: The equality of D((A + uB) k ) and D(A k ) will follow from the Kato-Rellich theorem 
( ||2T1 Theorem 1.4.2]). It suffices to check that for every k in N, u in R, e > 0, and every ip in 
D(A k ), there exists e < 1 and b £ such that 



IV. Tri-diagonal systems 




tend to zero. Then, for every k in N and u in 



',, that the sequence ( -f^- ) is bounded, 



\\((A + uB) k -A k )^\\ <e\\A k ^\\+b £ \m. 



(10) 
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Let us prove that B is bounded from D(A r+1 ) to D(A r ) for every integer r > 0. For every v 
in D(A r ), 

oo 

\\Bv\\ 2 r = J2^n\(Bv,<Pn)\ 2 
n=l 

oo 

n=l 

oo 

< E >ZKn-l I (0-1, W)| 2 + blMn, V)\ 2 + bl n+1 1 (0 n+1 , W)| 2 ) 



n=l 

00 / \ \ 2r h 2 h 2 

E « 2 -fF 1 ! (0n-l, «> | 2 + AT 2 ^f I «> | 2 + 



n=l 



v A„_i/ A n _ x A n 



/X \ 2r ft 2 
, \2r+2 / 71 \ u n,n+l\/i \|2 

+ K+l 7 "tJ IWn+uVjl . 

\A n+ i/ A n+1 

b 2 b 2 b 2 

Now for every e > 0, let n such that sup n > no ^f" 1 < sup n > no -^j 1 < e/3, and sup n >„ < 
|. Note that the sequence (A„) neN is non-decreasing. Then there exists C e such that 

no 

||^|| r 2 <EA 2r K^^n)i 2 + £ E \l r+ \^v)\ 2 <C £ \\v\\ 2 + e\\v\\ 2 r+1 . (11) 

n=l n>nrj — 1 

We prove (flOl) by induction on fc. For = 1 this is a consequence of (fTTI) with r = 0. The 
inductive step follows from the fact that 

(A + uB) k+1 - A k+1 = u{{A + uB) k B - A k B) + uA k B + ((A + uB) k - A k )A 

for every u in R and from inequality (fTTj) . ■ 

Proposition 10. Let (A, B) be a tri-diagonal system and let k be a positive integer. Assume that 
the sequence ( ^r 1 - ) is bounded, that the sequences ( "■ ra ~ 1 ) ( ^r 2 - ) ten( ^ to zero an d 

1 V An /nGN ^ V An / n gN V An / n gN 

the sequence |6 n ,n+i| ?^ 1 — 1 w bounded. Then (A, B) is k-weakly -coupled. 

Proof: For every ip in D(A), write ip = J2 < jLi%j<t>j where Xj = ((pj,ip). Since = 
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then 

oo 



ft U\A\ V, = & ( £ A/ / - r , / ', • . + A J. ; •' j • I hj.j ■ I •' ; 



oo 

!R ( jt^Xjfijbj+ijXj+i — Xjbj +1 jXj + i) + — \j)xj+ibj t j+iXj 



I I I "^3 "f" 1 1 



By hypothesis, there exists C such that |&jj+i|(A^ +1 — A??) < CXj for every j. Hence, |3ft(|A| fe ^, Bip)\ < 
CY!jLi tf\ x j\ 2 < C(|A|V; The equality of the domains follows by Proposition [9] ■ 

B. Estimates for tri-diagonal systems 

Lemma 11. Let (A, B) be a tri-diagonal system and I be an integer. Assume that the sequence 
( ^2+1 ) is bounded, that the sequences ( ) ( ^» ) ^ zero> an j 

V A n / nfN V A " / nfN V A " / neN 



neN 



ngN V A " /neN V A « / neN 

exists a positive integer k and < r < k/2 such that the sequences ^|6 n ,n+i| ^ r ^[ l — 1 

(iMv eN an< ^ (lvF"J eN are bounded. Then for every t > 0, /or every piecewise constant 
control u, 

ol 21+1 , t 

l\ j=l+1 \Jo 

where for j e N, L(j) = sup n m < i |6 r , 



J n,m 



Proof: Let K > 0. We prove the result for u piecewise constant of L 1 -norm smaller than 
K. For every e > by Proposition |4] there exists iV = N(K,e) > I such that ||T"(0i) — 
X^ N) {t, 0)0! || < e for every i > 0. 

Consider the solution if) : t i-> O)0i of ([|T/v} with initial condition 0i. Then ^(t) = 



e M(A °0i + Jo e(*- s ) AW M(s)S( 7V )^(s)c/s. Iterating / - 1 times we get 

m = e tA(N) 

i=1 Jo< Sj <-< Sl <t m=1 

+ / e (*--i)AW s W e («i-- 3 MW s (A0... e (-,-i--«)AW s (A0^ (aj ) A u ( Sm )d Sl ...d Sj Y 

7o< S; <-< Sl <t / 
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For the tri-diagonal structure of the system we have 
for every < Sj < ■ ■ ■ < s% < t and j < I — 1. Then 

(<f> l+1 , m) = e tAW [ (<f>i+i, e ^ A(N) B^e^-^ A(N) B^ ■ ■ . e ^-^ A(N) B^( Sl )) 

J0<si<—<si<t 

I 

Y[ u(s m )ds 1 ...dsi 

m=l 

Now, 

21+1 

sup \\B^e^- s ^ AW • • • c («-i)AW B^e^ A(N) cf> l+1 \\< 3 1 J[ L{j). (12) 
si,...,s ( e[o,t] i=«+i 

Then 



, 38+1 r 1 ,(tf\u(s)\ds) 1 2l + l 

m,^))i<3 ? n w) / n k^)i^i-.-^=3 a °' , ; n m 

j=l+l J0< Sl <-< Sl <t m=1 /! . i+1 



j=2+l - a ^' m=l i=(+l 

as a consequence 

l(0m,Tr(0i))|<3'— J] m+e, 

L j=l+l 

and the result follows as e tends to zero. ■ 
From a physical point of view, Lemma ITTI provides an estimation of the probability of energy 
transitions (in the spirit, for instance, of [18, Section X.12, Example 1]). 



Remark 6. In the case in which the diagonal of B is zero then equation (|12)) reads 

21+1 

sup \\BW e to-t-i)A™ B W . ..e^^B^e^-^^W < 2 l J[ L(j). 
si,...,s t e[o,t] j=i+i 

This gives the better estimate in this case | T^i) | < 2 l l[f=i+i W) ( Jo \ U ( T ) \ dr ) I 1 - 

C. Example: orientation of a rotating molecule in the plane II 

The system of Section IIII-CI provides also an example of tri-diagonal system. Recall that 
for this system, for every j,k in N, = k 2 , (0^,50^) ^ if and only if \j — k\ — 1 and 
(4>j, B(fij + i) = — i/2. We deduce a bound for the coupling constants from Proposition [TOl For 
every k in N, 

( X k \ 1 // 1 \ 2fc \ 2 2k — 1 

c k (A,B) < sup \{<p n ,B<p n+1 )\ -gl - l = sup - 1 + - -1= — — . 
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In particular c x {A,B) < 3/2 and, by ©, we obtain that || 7^17 (</>!) - XfafaOfrxfaW < e if 

/ I, || 3/2|M| rl \ 2 

A/y+i = (N + l) 2 > ^Mkif_ — j . 

The tri-diagonal structure allows to obtain better estimates on N. From Remark [6l we get 

(2K) N 2N+l K N 

\(<t>N + i,?n<t>i))\< { - 1 J- n m- " 



ATI 11 wy AT' 

iv - j=N+l iV - 

Therefore, by ©, for every e > 0, if iV is such that < 2£iV! then ||7T/vTJ< 

X^t, 0)0x11 <£. 

The second estimates is significantly better than the first one. For instance, if one has \\u\\ L i = 3 
and one desires e < 10~ 4 , the condition e(N+l) > \\u\\Lie 3 ^ 2 ^ Ll is false for every N < 2.7 10 6 
while the second condition, |H|^i +1 < 2eN\, is true for N = 14. 

D. Example: quantum harmonic oscillator 

The quantum harmonic oscillator is among the most important examples of quantum system 
(see, for instance, Il22~l Complement Gy])- Its controlled version has been extensively studied 
(see, for instance, 11231 . Il24l0 . In this example H = L 2 (R, C) and equation © becomes 

dib 1 

i-£(x, t) = - (-A + x 2 )tfj{x, t) + u{t)x^j{x, t). (13) 

An Hilbert basis of H made of eigenvectors of A is given by the sequence of the Hermite 
functions (</> n )neN» associated with the sequence (— iA n ) n£ N of eigenvalues where \ n = n — 1/2 
for every n in N. In the basis (0 n ) ne N 5 B admits a tri-diagonal structure: 



J 3i 



B<j) k ) 



—iy/k — 1 if j — k — 1 
—iVk if j — k + 1 
otherwise 

Proposition [9] and Proposition [lOl apply so that, for every k in N, the system (A, B) is /c-weakly- 
coupled and 

/ (n + l/2) fe \ 



Ck(A, B) < sup yfn 



(n - 1/2)* 



1 \ *4 / 1 



< supv^ 1 + r -l £ 1 + r 

3 fc+1 - 1 ^ 

< sup t 

2 n n — ~ 

< 3 fc+1 -l. 
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The quantum harmonic oscillator is not controllable (in any reasonable sense), however, 
the Galerkin approximations of (fT3l) of every order are controllable (see [23]). This is not a 
contradiction, indeed Proposition |4] states that there exists a size of the Galerkin approximation 
for which the trajectories of the infinite dimensional system can be approximately tracked by 
the Galerkin approximation, provided that there exists a uniform bound on the L 1 -norm of the 
control. As a matter of fact, there is no such bound for system (fT3l . 

To obtain an estimate of the order N of the Galerkin approximation whose dynamics remains 
e close to the one of the infinite dimensional system when using control with L 1 -norm K, one 
could use Proposition 0] with k = 2, r = 1, d = 1, and ||0i||i = 1/2. The resulting bound, as 
given by d8]), 



is however very weak. Like in the example of Section IIV-CI the tri-diagonal structure of B 
allows better estimates. Using Remark [6l we find that \\X^ N \t,0)(j)i — 7TArT"0i|| < e provided 
IMU 1 < K and 



For instance, if K = 3 and e = 10~ 4 , this is true for iV = 420, while COO) is false for N < 10 



In our study we focused on the notion of weak coupling. We have established some inter- 
esting consequences in control theory and numerical simulation which applies to common and 
interesting physical models. 

However, our assumptions are not optimal and in forthcoming works we expect generalization 
to rough control such as Dirac impulses. We hope that we will include systems with continuous 
spectrum in the scope of such technology. 



Proof of Lemma [7} Without loss of generality we can assume that the operators \A\ 
and \A'\ are positive and invertible. Let (0 n )neN and (0^) neN be unitary basis of H made 




V. Conclusion 



Appendix 



Proof of Lemma CD 
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of eigenvectors of A and A' respectively. Then A„0 n = \A\<j) n for n G N and Z)(|A| S ) = {ip G 
H : E ieN Af |(^,^)| 2 < +oo}. Similarly, we can define A^ and D(\A'\ S ). 

Since _D(|A| fc ) C D(\A'\ k ) and by the closed graph theorem, we deduce the existence of 
C k > such that for every ip G D(\A\ k ) 



E A « m€)\ 2 <C k Y,\^n\ 2k m<Pn) 



so that 



,2k 



E A n 



Y,{^Aj){<t>jA'n) 



For all V G let ^ in # such that ip = \A\- k ip = E Aj*^, Then, for all ip £ H, 

we have 

,2fc ■ 



E A n E 0/) (<f>l, <P'n) E A i ^ #i) < C<A 



(15) 



n ; j 

and the equality holds for /c = and C = 1. Consider ip e H and 

: * = » + iy E At +j2i/ (l^r s+i ¥, 14—*$ 

n 

where, for every 2 in C, = Ej A|(^, 4>j)4>j- Then, by (fT?T) for s = and s = k we have 



/:*(* + iy) <C S |||A| 



iai— ^il < a 



If ^ is finite linear combination of the vectors {0j}jeN then the function /- analytic on the strip 
{z G C : < tflz < k} and continuous on its closure as uniform limits of a partial sum on n. 
Since it is bounded on the boundary, by Hadamard three-lines theorem [fT8l Appendix IX.4], it 
is bounded on the strip, and, moreover, log(sup R2=s is a convex function of s G [0, k]. 

So that for s G (0, k), we obtain 



E A n KV>, <0 2S < Cfc E \ X n\ 2S \(lp 



1 rn/ I 5 



and by density _D(|y4| s ) C D(\A'\ S ). The hypothesis and the proof being symmetric in A and A' 
we have actually the equality. ■ 
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Appendix 
Proof of Proposition [2] 

In this Appendix, we use regularization techniques to provide a proof of Proposition [2] 
Proof of Proposition^ Note that for every u E R p , D(\A + EiUiB^/ 2 ) = D(\A\ k l 2 ), 
the function \A\ k / 2 e^ A+ ^ UlBl) i) is in C(R,H) and for every e > the function \A\ k / 2 (e(A + 
EiUiBi) + l^V^+E^Vo is in C\R,H) whenever ip E D(\A\ k l 2 ). 

If t i— )■ ^>(t) is the solution of © with initial condition ip in Z}(|A| fc//2 ), the real mapping 
/ : t (->• (|A| fc, 0(*)) ^(^)) is absolutely continuous from R to R. We make a regularization 
to obtain extra regularity, we introduce /J : t + uj" 1 Bi) + l)' 1 ^^), (e(A + 

Y^i^i^Bi) + From the functional calculus ll25l Theorem VIII.5] the sequence fi is 

pointwise convergent to / as e tends to 0. 

The function /J is absolutely continuous from R to R and it is differentiable on the interval 
(tj-i,tj), for every t E (tj-i,tj), 

j t fm = j t (\A\\e{A + Y.< lB i) + 1)"VW, (^ + Ert) + 1)"V(*)> 

= <|A|*((i4 + E + (A + t u i(t) B i) HA + E ut'Bt) + 1)"V(«)> 

z Z=l z 

+ (l^l fc (^ + EW" 1 ^) + i)- 1 ^ + E uMBMt), «a + EW _1 A) + i)"V(*)> 

Z Z=l z 

= 23R(|A| fc (e(A + E < l Bi) + 1)"V(*), (A + E «i(*)5i)(e(i4 + E + 1)~V(*)> 

z 1=1 i 

= 2 E «i(f)»<l A\ k {e{A + E + 1)"V(*), + E + I) -1 ), 

Z=l z z 

and since (A, Bi, . . . , B p ) is A;-weakly-coupled, 

|^/|(t)|<2c fc (A,S 1 ,..., J B p )x 



x E hz(t)||(|A| fc (e(A + E + l)-ty(t), (e(i4 + EW^) + 1)"VW>I 
Z=l z z 

<2c fc (A, J B 1 ,..., J B p )E| M/ (t)l/ e J W- 



z=i 



Gronwall's lemma implies that = {\A\ k (e(A + Y,i uf^i) + (e(A + £,u{~ 1 fl,) + 

l) -1 ^(*))<e '^mVi 1 111 ' fi(tj-i). Passing to the limit e to 0, this gives f(t) = 
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(\A\ k if)(t),ip(t)) < e 2cfc ^ A ' Bl ''"' Bp ^i=ii'%-i' U! '^ d ' r j^ j ._ 1 ^ ^n immediate iteration concludes the 
proof. ■ 
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